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Abstract 

We  analyze  the  properties  of  electromagnetic  metamaterials  with  anisotropic  constitutive  parame¬ 
ters.  Particularly,  we  analyze  the  so-called  Radial  Photonic  Crystals,  which  are  radially  periodic 
structures  verifying  the  Bloch  theorem.  This  type  of  crystals  can  be  designed  and  implemented  in 
acoustics  as  well  as  in  electromagnetism  by  using  anisotropic  metamaterials.  They  were  originally 
proposed  in  acoustics,  but  similar  functionalities  are  here  proven  in  the  electromagnetic  domain 
together,  for  the  first  time,  with  an  analysis  of  their  behavior  operating  in  the  microwave  regime. 
Finally,  we  present  a  complete  discussion  concerning  their  electromagnetic  properties  together 
with  a  practical  design  of  such  a  microwave  structure. 

1.  Introduction 

Anisotropic  phenomena  have  been  largely  investigated  both  in  acoustics  and  electromagnetism.  Con¬ 
trol  of  anisotropic  propagation  is  a  key  condition  to  allow  a  number  of  applications,  which  can  cover  a 
broad  spectrum  of  areas  from  directive  antennas  to  cloaking  devices.  In  this  context,  periodic,  semi- 
periodic  and/or  multilayered  microstructures,  depending  also  on  the  desired  spatial  configuration  (ID, 
2D  or  3D),  have  been  studied  to  satisfy  different  types  of  targets.  These  studies  are  usually  based  on 
the  analysis  of  the  desired  propagation  characteristics  at  the  macroscopic  level,  and  of  the  elementary 
constituent  cells  at  the  microscopic  level.  In  particular,  the  selection  of  arrangement  and  unitary  cell 
element  geometry  is  inherently  linked  to  the  targeted  spatial  configuration,  wave  polarization,  opera¬ 
tion  frequency  and  dispersive  characteristics,  and  isotropic  or  anisotropic  behavior,  among  others.  We 
analyze  and  compare  the  wave  propagation  characteristics  of  the  so-called  Radial  Wave  Crystals 
(RWCs),  [1],  or  Radial  Photonic  Crystals  (RPCs)  in  the  electromagnetics  field.  This  type  of  metamate- 
rial-inspired  microstructure  has  the  original  characteristic  of  being  invariant  under  translation,  by  im¬ 
plementing  radially  dependent  constitutive  parameter  functions  in  multilayered  systems.  Although  the¬ 
se  structures  were  first  considered  for  acoustic  waves,  clear  analogies  can  be  established  between  both 
application  fields,  [2].  For  a  2D  configuration,  the  fundamental  characteristic  of  these  microstructures 
derives  from  the  periodicity  of  the  radial  profile  of  the  anisotropic  constitutive  parameters.  In  this 
work,  this  is  first  analyzed  by  means  of  the  proper  design  equations,  together  with  the  dispersion  dia¬ 
gram  information  and  related  propagation  parameters.  Then,  finite  size  structures  based  on  these  con¬ 
cepts  are  studied,  the  excitation  of  RPC  structures  with  line  sources  is  assessed  and  some  potential  ap¬ 
plications  are  described.  Finally,  a  practical  design  of  a  sample  microwave  microstructure  is  per¬ 
formed  taking  into  account  some  simplifications  that  allow  the  feasibility  of  the  device.  These  limita¬ 
tions  are  imposed  on  the  number  of  constitutive  parameters  and  layers  implemented.  This  paper  ends 
with  a  conclusion  section  summarizing  the  main  findings  of  this  work. 
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2.  Radial  Wave  Crystals  for  electromagnetic  waves 


Starting  from  the  scalar  Helmholtz  wave  propagation  equation,  for  the  radial  part  and  anisotropic  pa¬ 
rameters  and  TMZ  waves  (Ez  field),  we  have: 
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where  sz(r %  ji^r)  and  jur(r)  are  radially  dependent  constitutive  parameter  isotropic  functions  and  q  is  a 
constant.  Equation  (1)  is  a  Bessel  equation  of  order  q  that  admits  solutions  invariant  under  translations 
of  type  r  — ►  r  +  nd,  [3].  A  further  analysis  of  the  problem  allows  calculating  a  dispersion  relation  for  a 
multilayered  structure  made  of  alternating  layers  of  types  a  and  b ,  with  periodicity  d  =  da  +  db  as: 
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Figure  1(a)  displays  the  radial  dispersion  diagram  for  a  RPC  with  constitutive  parameters  verifying  the 
previous  conditions.  Figure  1(b)  shows  a  set  of  constitutive  parameters  allowing  the  implementation  of 
a  ‘reduced’  parameter  RPC  structure  with  4  layers  (only  sz(r)  and  jud(r)  are  implemented  and  jur(r)  =  1). 
This  reduced  profile  set  behaves  and  is  equivalent  to  the  one  of  the  ideal  structure  of  Fig.  1(a). 


Fig.  1:  (a)  Dispersion  diagram  for  the  first  5  modes  of  a  Radial  Wave  Crystal.  The  RPC  is  made  from  a  multi¬ 
layered  structure  of  alternate  layers  of  types  a  and  b  with  constitutive  parameters  (//ra’7,  ju0tl,  £za1)  =  (0.347r/d, 
0.08 r/d,  0.143r/d)  and  ( ju rb\  fdeb,  s^'1)  =  (0.5 r/d,  0.04 r/d,  0.1  r/d)  and  periodicity  parameter  d  =  da  +  db.  Field 
patterns  correspond  to  the  design  in  the  right  panel;  (b)  Reduced  model  parameters  for  a  four  layers  practical 
implementation  as  a  function  of  the  radial  distance  (center  of  the  RPC  located  at  r  =  0  mm).  Angular  permeabil¬ 
ity  (solid  line)  ju^r)  follows  a  stair-like  profile  bounded  by  layers  a  and  b  analytical  parameter  functions  (dotted 
lines);  permittivity  (dashed  line)  sz  =  3.4  for  all  the  RPC  shell;  radial  permeability  jur  (not  plotted)  is  constant  and 
equal  to  that  of  the  background  (jur  =  1). 


3.  Practical  realization  of  a  Radial  Photonic  Crystal 

For  a  practical  implementation  of  the  constitutive  parameters  we  have  selected  a  unitary  cell  composed 
of  a  resonant  split  ring  resonator  (SRR)  particle.  This  unitary  component  has  reduced  bianisotropy  ef¬ 
fects  compared  for  example  with  other  SRRs.  This  is  important  due  to  the  fact  that  it  is  not  desired  that 
the  array  of  resonators  implementing  the  /u^r)  function  interferes  with  the  /ur(r)  function  (which  in  our 
case  should  be  neutral  jUr(r)  =  1).  Generally,  the  permeability  of  an  array  of  SRRs  follows  a  Lorentz- 
like  model  with  the  resonant  frequency  separating  positive  and  negative  values  of  effective  permeabil¬ 
ity.  With  the  proper  design  of  the  geometric  dimensions  of  the  SRRs  it  is  possible  to  define  the  perme¬ 
ability  response  at  a  specific  frequency  (the  design  frequency).  An  optimization  process  is  required  to 
match  the  desired  permeability  and  permittivity  values  at  the  desired  operation  frequency  (3.8  GHz,  q 
=3  resonance  in  our  case).  Figure  2  displays  a  fabricated  device  based  on  parameters  of  Figure  1(b), 
together  with  the  unit  cell  and  the  relevant  geometric  dimensions.  We  intend  to  verify  the  interaction 
of  a  line  source  in  close  location  to  the  device. 
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Fig.  2:  (a)  Fabricated  4  layers  prototype  made  of  Split  Ring  Resonators  and  a  Rohacell  (sr  ~  1)  foam  support,  (b) 
Unit  cell  configuration  to  implement  the  angular  permeability  profile  with  characteristic  dimensions. 

The  real  part  of  the  E-field  patterns  in  the  region  around  the  RPC  are  displayed  in  Fig.  3.  The  RPC  is 
centered  on  the  origin  of  the  radial  coordinate  system  and  the  line  source  is  located  at  position  x0  =  yo 
-  95  mm.  Figure  3(a)  shows  the  result  for  an  analytical  model  simulation  (Comsol  software).  Figure 
3(b)  shows  the  simulation  results  for  the  SRR  implemented  RPC  (cut  plane  plot  of  3D  simulation  in 
Ansoft  HFSS  software),  while  Figure  3(c)  shows  the  measured  response  of  the  fabricated  device.  Col¬ 
or  scale  represents  real  part  of  E-field  with  same  range  in  all  plots.  From  this  comparison,  it  is  clear 
that  the  q  =  3  mode  is  excited  in  the  RPC  shell  at  the  predicted  resonance  frequency,  f  =  3.80  GHz. 
Field  maxima  and  minima  are  closely  comparable  in  both  simulations  and  the  measurement.  Orienta¬ 
tion  of  the  mode  pattern  towards  the  source  is  observed.  Within  the  RPC  shell,  overall  the  field  pat¬ 
terns  closely  follow  the  ‘envelope’  corresponding  to  the  q  -  3  field  pattern. 


(a) 

(b) 

(c) 

r<3> 

A 

'*w 

%«> 

*  %  F 

» 

Fig.  3:  Real  part  of  the  E-field  of  a  line  source  externally  illuminating  a  RPC  shell  at  3.8  GHz  (a)  analytical  defi¬ 
nition  of  the  RPC  constitutive  parameters  (Comsol  simulation);  (b)  ab  initio  simulation  of  the  RPC  as  in  Fig.  2(a) 
(HFSS  simulation);  (c)  measured  results  on  the  fabricated  device. 

4.  Conclusion 

Radial  photonic  crystals  are  ‘very  ordered’  microstructures  that  have  been  analyzed  by  means  of  their 
electromagnetic  constitutive  parameters  profiles  showing  their  possible  application  as  resonant  struc¬ 
tures.  Their  combination  with  radiation  sources  allows  the  analysis  of  the  interactions  and  points  out 
possible  applications  for  example  as  beam-forming  shells.  Their  application  as  frequency  and  position 
sensors  can  also  be  anticipated.  The  behavior  of  these  microstructures  was  demonstrated  with  a  proto¬ 
type  device  operating  at  3.8  GHz,  whose  response  showed  excellent  agreement  with  our  model. 
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